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The dynamics of random weakly nonlinear waves is studied in the framework of vibrating thin 
elastic plates. Although it has been previously predicted that no stationary inverse cascade of con¬ 
stant wave action flux could exist in the framework of wave turbulence for elastic plates, we present 
substantial evidence of the existence of a time dependent inverse cascade, opening up the possibility 
of self organization for a larger class of systems. This inverse cascade transports the spectral den¬ 
sity of the amplitude of the waves from short up to large scales, increasing the distribution of long 
waves despite the short wave fluctuations. This dynamics appears to be self-similar and possesses a 
power law behavior in the short wavelength limit which is significantly different from the exponent 
obtained via a Kolmogorov dimensional analysis argument. Finally, we show explicitly a tendency 
to build a long wave coherent structure in finite time. 

PACS numbers: 05.45.-a, 47.27.E-, 46.40.-f, 62.30.-|-d 


I. INTRODUCTION 

Oscillating random waves are present in a myriad of situations in nature, displaying a large variety of scales and 
exhibiting turbulent-like behavior, the so-called wave turbulence [H-Q. Of particular interest are the oscillations 
over the surface of the sea, Rossby waves in atmospherical science, nonlinear optics, the plasma oscillations and the 
vibration of elastic bodies such as piano strings, timbals, or more complex singing bowls, bells or gongs. Because of the 
intrinsic nonlinearity of the basic underlying physics of these problems, and because of the randomness of the phase of 
the oscillations, only a statistical description seems reasonable. The weak turbulence theory provides such a statistical 
description for the asymptotic long time behavior of the spectral wave amplitude, in the case where nonlinearities are 
small. In particular, it describes the energy transfer among the different modes in agreement with the conservation of 
the total energy of the waves. More precisely, this wave turbulence theory provides kinetic equations for the long-time 
evolution of the spectral amplitude for dispersive wave systems [H-Q. In the present context of small nonlinearities, 
we will use interchangeably wave turbulence theory and weak turbulence theory and refer to it as WTT. Remarkably, 
such kinetic equations exhibit stationary solutions corresponding to equipartition or constant flux cascades of the 
energy, namely the Rayleigh-Jeans solution and the Kolmogorov-Zakharov (KZ) spectrum respectively. The search 
of the KZ spectra has motivated exhaustive studies these last fifty years Q, regaining recent interest, because of the 
parallel development of new theoretical and experimental findings. Among them, we mention the cases of surface 
capillary waves in, surface gravity waves [1, 0 and elastic waves of thin plates [M3 for instance. While such 
dynamics corresponds usually to a direct cascade of energy towards the small scales, the formation of large scale 
structure can sometimes be observed. This is the case in particular for Bose-Eintein condensation or in two 

dimensional hydrodynamic turbulence [l5j| , where the self organization process is a consequence of an inverse cascade. 
This inverse cascade transfers some quantity (particles, enstrophy, wave action) from the small scales toward the large 
scales leading to the formation of coherent structures. The formation of an inverse cascade in different systems has 
always been related to the existence of a conserved quantity at least at the weakly non linear level. For instance, 
when using the Gross-Pitaevskii or nonlinear Schrddinger equation to model the Bose-Einstein condensates, WTT 
predicts the existence of an inverse cascade of mass (a conserved quantity). Similarly, in the case of surface gravity 
waves, an inverse cascade of the conserved wave action is deduced and numerically observed da El- Finally, we want 
to emphasize that besides the description of the stationary solutions of the dynamics, the kinetic equations that are 
deduced by the WTT give a very good framework to investigate non stationary situations involved in wave systems 
such as transitory or decaying regimes that often lead to self-similar dynamics [l9l - [^ . The goal of this paper is to 
show, using the elastic vibrating plate, that the formation of an inverse cascade does not generally require a conserved 
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quantity, opening up the possibility of self organization for a larger class of systems. In particular, nothing prevents 
the existence of a time dependent inverse cascade that would transfer wave action from short scale to large scale. 

The paper is organized as follows: Section|lll introduces the dynamical version of the Fdppl-von Karman equations 
which provide the basic nonlinear equations for vibrating elastic plates, containing both bending and stretching. 
Then we summarize the main findings of the WTT of a vibrating plate, in particular the concept of the Kolmogorov- 
Zakharov spectra. Section IIIII presents the numerical simulations of a vibrating plate which is forced only at very 
short wavelength, displaying a striking inverse cascade of wave action. Section IIVI analyzes the numerical evidence of 
a self similar evolution which suggests a blows-up in finite time. Finally, we conclude with an overall discussion of the 
problem. 


II. WAVE TURBULENCE THEORY OF VIBRATING ELASTIC PLATES 
A. The Foppl—von Karman equations for elastic plates 

Vibrating elastic plates offer, perhaps, the most suitable weakly non linear wave system. It is studied within the 
framework of the dynamical version of the Foppl-von Karman equations which model the dynamics of the out of 
plane displacements of the plate. We shall use the same notations as in Ref. @ , but we write them in dimensionless 
units. We choose I = h/-^3(1 — cr^) as the unit of length, and as the unit of time. Here h is the thickness of 
the elastic sheet, the material has a mass density p, a Young modulus E and its Poisson ratio is a. In these units the 
equations read: 


^2 A 1 

^ = --Arctic, x}; w 

A^X = -^{C,C}- (2) 

The out of plane displacement of the plate in physical units is thus and the Airy stress function is 

Equation (I2|) for the Airy stress function x(a;,y,t) may be seen as the compatibility equation for 
the in-plane stress tensor which follows the dynamics at the lowest order [^. The characteristic size of the plate is 
L, thus the dynamics of a free plate is governed by a single dimensionless parameter: A = j = -^3(1 — cr^)^, which 
is typically of the order of 10^ up to 10^. A = dxx + dyy is the usual Laplacian and the bracket {•, •} is defined by 
{/, g} = fxxPyy + fyydxx — ‘i-fxydxy, which is an exact divergence, so that equation ([TJ preserves the momentum of 
the center of mass, hence, du f C(x,y,t) dx dy = 0. Moreover, the total energy: 

^ = / + ^(ACf-^(Axf-^x{C,C}) dxdy, (3) 

is also conserved by the dynamics (HED- Finally, small plane waves perturbations (( ~ x ^ = {x,y)) of 

a plane plate are dispersive with the usual ballistic behavior of bending waves, that is Wfc = [^ . 


B. Wave turbulence equations for the spectral densities 

As already discussed in Ref. Q, equations (I1I2I) exhibit a Hamiltonian structure which is easy revealed in Fourier 
space, defined by Ck{t) = ^/ C{x,t)e^'^'^d‘^x, with = C-fc- The Hamiltonian structure allows one to performs a 
canonical transformation 

Ck = —7=={A-k + A-k) ( 4 ) 

V2ujk 

which lead to a diagonalized form of the wave equation: 

dAk 

—— -I- iuJkAk = iNaiAk), (5) 

dt 

where A^ 3 (-) abbreviates the cubic nonlinear interaction term given explicitly in Ref. Q. 

The WTT describes the long time statistical behavior of weakly non linear random waves. The analysis is based 
on an infinite hierarchy of integro-differential equations for the cumulants of the canonical variables which maybe 
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deduced directly from ([5]). In the weak wave amplitude limit, a multiscale asymptotic expansion of these hierarchy of 
equations provides a rational scheme for solving every cumulants EH!- As a result, the second order cumulant 

+ ^ 2 ), ( 6 ) 

is shown to control the long-time dynamics of the wave system, where nk is called the spectrum of the wave. Other 
second order cumulants vanishes in the weak amplitude (long time) limit , in particular 

0 & ^ 0 . 

The asymptotic perturbation scheme of this theory provides at first order a nonlinear frequency shift to the linear 
waves, leading to an effective oscillation frequency = cok + . This correction due to weak nonlinear 

effects is a function of the mean spectral density nk{t) ([SI) and it reads (in the dimensionless units) [ 1 ^ : 






riq qdq - 



-^n„ qdq 
Q 


(7) 


Notice that this frequency correction was also obtained by considering a limited number of nonlinear interactions [^ . 
In addition, in the WWT this frequency shift © is useful to quantify the nonlinear effects, and by consequence the 

validity of the WTT. Indeed, the ratio uj^'^ /cuk indicates the relative importance of the nonlinear term with respect to 
the linear behavior. The uniform validity of the WTT requires that this ratio should satisfy 1 for all the 

wave numbers. If this number is of the order of unity, wave turbulence is no longer valid, at least for the concerned 
scales. 

At the next order, WTT provides a kinetic equation that governs the mean spectral density evolution nk(t), which 
reads Q: 


dt 


Hk — — f‘27r J‘ dlCidk2dks \J—k^ki;k2^k^\ 

X nk^Uk^Uk^nk { — - — - — - —] - ki - k2 - k3)S{uJk - Siojki - S2UJk2 - ssOJks)- (8) 

\nk nk, nk2 UkJ 


The coefficient J-k,kiM,k 2 c comes from the fourth order nonlinearities in the total energy ®, and they 

are given explicitly in Ref. [ 8 |. The details of this scattering function J is not needed here and we omit to write it for 
the sake of simplicity, but, for the purpose of this work, we only have to notice that it has zero degree of homogeneity 
in fc, that is 


d~Xk,Xki',Xk2,Xk3 — J—kMi\k2.k3- 

Finally, we mention that this function J-fc.fc 1 .fc 2 .fc 3 vanishes as fc —^ 0, so that the spectrum, does not vary at fc = 0, 
in agreement with the original plate equations © and ([2|). If nk=o = 0 at t = 0, then the spectrum vanishes at k = Q 
for all time. 

The spectral dynamic described by the kinetic equation ® corresponds to four-waves resonances which enforce the 
energy and momentum conservations in each interaction. However, in contrast with the case of diluted gases where 
each collision preserves the number of particles, the total number of waves involved in the interaction is not formally 
conserved for vibrating plates. Moreover, it is interesting to notice that while the total energy ® is conserved by 
the original FoppNvon Karman equations ©3), the kinetic equation preserves only the quadratic part of the total 
energy, namely, A 2 = / (^(^tC)^ + |(AC)^) dxdy. 

The sum in m rules for Si = ±1, that is, the kinetic equation ([5|) contains 8 terms. Among them, the one 
corresponding to all Si equal is not resonant and it thus vanishes. Three other terms (all identical by symmetries) 
correspond to interactions of two waves coming-in and two waves going-out, so that the total number of waves is 
preserved by this interaction. We will refer to these terms further on as the 2 0 2 resonant case. Finally, there are 
four other interaction terms corresponding to one(three) wave(s) coming-in and three(one) waves going-out, which 
are referred as the Sol resonant case. In this later case, the total number of waves is not preserved formally by 
the four waves interaction. We denote these two different interaction terms by C 22 and C 13 , respectively, so that 
C[nk] = C 22 [nk] + Cisluk]- This non wave action conservation represents a major difference with most of the known 
four-wave interaction systems such as surface gravity waves fl^ . or nonlinear optics [ 2 ^. To our knowledge the only 
known physical systems that exhibit these two kind of interactions (202 and Sol) are the symmetric capillary 
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waves at the interface between two fluids @ and the elastic plates Q. 

In conclusion, although the kinetic energy £ = J LOknk{t)cPk, is preserved by the dynamics ([5]), the wave action, 
A/” = / nk{t)(Pk, is not. 


C. Kolmogorov-Zakharov spectra 

Although the wave action is not preserved by the dynamics, local conservation equations can be deduced from the 
kinetic equation. Indeed, the change in time of the energy spectral density E{k) = 27rfcwfenfe, can be written, after 
(0), as 


d d 

-—Elk) = ——P{k) where P{k) = 2Tr / u;qC[nq]qdq (9) 

(it drZ J 

is the energy flux, which depends, in principle, explicitly on the wavenumber k and t. Similarly, the wave action flux 
Q{k) maybe defined via the wave action spectral density: N(k) = 27rfcnfc, through: 




with 


Q{k) = 2 tt [ C[nq\qdq. 
Jo 


( 10 ) 


As the energy flux, the wave action flux may depend on the wave number and time. Notice that the special writing 
of equations m may induce the wrong impression that J N(k)dk is conserved by the dynamics, but this is not so, 
because <5(0) ^ Q{oo). Wave turbulence theory predicts a class of exact power-law solutions of the kinetic equation 
(0), found by Zakharov [I|, which keep the fluxes constant. More precisely a (direct) energy cascade is found for 
which the energy flux P is constant. Similarly, if the wave action is conserved by the dynamics, then, an (inverse) 
wave action cascade corresponding to a constant wave action flux Q can be exhibited. These solutions are named the 
Kolmogorov-Zakharov (KZ) spectra, because Zakharov’s findings are in close relation with the Kolmogorov scaling 
arguments used in fluid turbulence. 

Taking an arbitrary power law solution, = A k~'^^, introducing this into the collisional operator C[nk\ one readily 
gets: 


C[Ak-'^^]= A^I{x)k^-^^. 

Here is a pure function, which depends only on the exponent x, and whose expression has been explicitly written 
in Ref. Q. Following (0) and (fT0|) . the fluxes are then given by: 

P = A^ ^^(^) fc6(i-^) and Q = (II) 

6 (x-l) ^ (2-3x) ^ ’ 

Constant energy or wave action fluxes are obtained if the exponents take the values x = Iorx = 2/3 respectively. 
Such arguments guarantee only that the scaling of the solution is consistent with the collisional operator. However, 
since the denominator vanish for those exponents, one needs in addition that the collisional operator vanishes in order 
to obtain stationary solutions of the kinetic equation. This condition is not satisfied for elastic plates, where the 
inverse cascade of wave action, x = 2/3 is not a root of /(x). Only the terms corresponding the 2 0 2 resonances 
vanish for x = 2/3, while the terms due to 3 O I resonances do not. Then, the flux Q formally diverges and the KZ 
solution is not valid. 

On the other hand, /(x) has a double root for x = I due to a special degeneracy (/(x) vanishes quadratically near 
X = 1, 1{x) ~ (1 —x)^). It indicates that both the Rayleigh-Jeans and the KZ solutions exist for x = 1. In practice, the 
resulting flux is zero, thus a logarithmic correction should be included on the final spectrum 0. Therefore, the above 
considerations imply that only one cascade is guaranteed, namely a direct cascade of energy toward the small scales, 
which remarkably is not a simple power law. This KZ spectrum, predicted by [^, reads (the numerical pre-factor is 
discussed in Ref. El): 


^Mrect ni/3 ln^^^(fc^/fc) 
''■fc ^ ^2 


( 12 ) 


where P is the energy flux, and k^ a cutoff scale. These solutions have been observed in numerical simulations 
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performed with an ad-hoc dissipation concentrated at small scales only However, experimental observa¬ 

tions d, [^.[^ [2^ present a slightly different behavior for the direct energy cascade, which is understood as follows. 
In Ref. [30], it is shown that the dominant low frequency dissipation rate of the damping suppresses the existence 
of a window of transparency in the wave number range probed by the experiment. It turns out, that the stationary 
spectrum for a vibrating plate comes from the balance between the kinetic collision integral, the forcing and d amp ing, 
displaying no simple KZ scaling spectrum, and exhibiting a strong dependence on the damping mechanisms [30|. 

We conclude this section with the following remark. If one neglects the 3 1 resonances in the kinetic equation, 

that is if one imposes Ci3[nfc] ~ 0, then an inverse cascade of wave action with constant flux can exist in addition to 
the direct cascade of energy. This inverse cascade of wave action (from small scales to large scales) reads: 


^inverse 




1 

^4/3 ’ 


(13) 


where Q < 00 is identified as the (constant) wave action flux. This inverse cascade transfers wave action between 
modes, a self-organization process which may lead to the formation of coherent structures and eventually to the 
breakdown of the WTT [3, [2^ [3l| . 

The goal of this paper is therefore to investigate the plate dynamics by forcing the vibrations at small scales only in 
order to observe the genuine transfer of wave action towards large scales despite the presence of the 3 e-)- 1 interactions. 


III. MANIFESTATION OF AN INVERSE CASCADE 

We solve numerically the coupled set of dynamical equations (HHH) using a pseudo-spectral method which takes 
advantage of the linear wave dynamics in Fourier space. Formally, the equations dim read, in Fourier space: 

Cfe = — WfcCfe + ./VLfe — HfcCfc -I- /fc, (14) 

where NLk stands for the Fourier transform of the nonlinear term including equation Q, Dfe represents a linear 
damping and Ik the forcing in spectral space. Finally, the temporal integration is performed in the Fourier space 
using a second order Adams-Bashford scheme. Standard dealiaising technique for cubic nonlinearity has been tested 
in previous works @ with no qualitative changes in the results so that no dealiaising is used in the present simulations. 

In the present work we use periodic boundary conditions, which are the natural framework to investigate the features 
of the wave turbulence, the number of modes ranging from 512^ up to 2048^, with a mesh size dx = 1/2, leading 
to the spectral ultraviolet cut-off kc = tt /dx = 2'k. For numerical stability, the time step used for the simulations is 
dt = 0.02 unit time. To observe the dynamics towards large scale, we force and dissipate the system at small scales 
only. The dissipation is given by Dk = —rj{k^ — k‘^)H{k — kd) where rj is the amplitude of the damping, H{-) is the 
Heaviside function and 2tt/ kd the characteristic scale below which the only dissipation acts. The forcing will be non 
zero only in the finite range [ki — 6i,ki + 5i\ where 

4= A,--e 

Here ki is the characteristic scale of the forcing, Si and Ai its the width and amplitude respectively. The angular 
variable 0k(t) is a random phase taken in the interval [0, 27r], which also changes randomly in time. Notice that this 
process injects both energy and wave action around ki since one cannot separate them formally. To illustrate this 
inverse transfer mechanism, we take kd = ki + Si so that the inertial range for an energy cascade vanishes and the 
low frequency (large scale) inverse transparency window is the largest possible available. Finally, we have checked 
numerically that the results do not depend on the details of the dissipation at small scales. 

Although it is required to dissipate the energy at small scales (large k) to reach numerically a stationary state, we 
have realized distinct numerical simulations with a sink and without it located near fc = 0, and we conclude that it is 
not required to absorb nor dissipate the energy (nor the wave action) at the large scale (small k) to reach numerically a 
stationary state in the time scale of the simulations. This can be explained firstly by the fact that energy is eventually 
dissipated at small scales leading to the general balance between the injected and the dissipated energy. For the 
wave action, since this quantity is not conserved by the dynamics, everything works as if the wave action is formally 
absorbed by a sink at fc = 0 (which is a neutral mode) so that there is no need to add such an absorption term near 
fc = 0 in the dynamics. 

Fig. [1] shows the snapshots of the plate deformation at six distinct times of the evolution. Notice the apparent 
formation of a coherent structure which at the end oscillates at the largest possible mode. This coherent structure 
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FIG. 1: Snapshots the out of plate displacement (^(x,y,t) obtained by the numerical simulation of the Foppl-von Karman 
equations (I1I2II at six different times : a) t = 22 x 10^; b) t = 50 x 10®; c) t = 62 x 10®; d) t = 82 x 10®; e) t = 162 x 10®; and 
t) t — 202 X 10®. The injection is made at small scales, with ki = 4.5 and Si = 0.5. Numerical dissipation acts at smaller scale, 
starting at the end of the injection range {kd = 5). The amplitude of the injection is Ai = 0.0001. The system size is 1024^ 
units with 2048^ modes and kc = 27r. Note that the vertical scales differ from one figure to the other, growing from a) to f). 
Despite this change of scale, the amplitude of the small spatial scales can always be observed. Notice the formation of a large 
scale structure as time increases. 


appears as a consequence of the long time evolution which is mostly characterized by the largest modes of oscillation 
of a plate plus small fluctuations. 



FIG. 2: a) Evolution of the energy density E/L?", and, b) the wave action density J\f/L'^ with time. The parameters are the 
same than those of Fig. [T] The time of the different snapshots of Fig. [T]are indicated on the curves. 

Fig. [2] shows the numerical evolution of the energy and the wave action with time for this numerical simulation. 
After a transitory regime where both quantities vary, we observe that a quasi-stationary regime is reached above 10® 
unit time approximately. 

This dynamics can also be inspected within the evolution of the wave spectrum which is defined following eq. o 
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FIG. 3: a) Angular average of the spectra, n*,, as a function of wavenumber k in log-log scale, at different time steps starting at 
f = 35 X 10®. Subsequent spectra are labeled according to the snapshots of Fig. [T] The straight line indicates a power law 
as a reference guide, b) Under the same conditions log-log plot of the wave action flux Q{k,t) as a function of wavenumber for 
the same times. 


and dni by: 


nk{t) = UJkilCkf). 

Fig. [3]-a shows the averaged value of the wave spectrum over the angle in the fc-space.The forcing creates a wave 
action flux towards k = 0 that “fills” the spectrum at large scale. This can be observed on Fig. [3]-b which present 
the wave action flux m at different time, computed explicitly as a sum over discrete modes Q{k) = kdtnu- 

(Note that by definition <5(0) = 0.) For large time (again, above 10® time units) the spectrum tends asymptotically 
to a stationary form that exhibits a power law with an exponent surprisingly close to the hypothetical 4/3-inverse 
cascade exponent m which is forbidden by the 3 O 1 interactions. Similarly, the wave action flux Q converge to an 
almost constant value. Notice however, that the dynamics is not steady but only in a “quasi”-stationary regime since 
the large scale modes still exhibit a slow dynamics. 

In the following, we will show that the evolution of this amplitude spectra can be decomposed in two distinct 
regimes in time, both being dominated by weakly random waves. The first regime displays a self-similar behavior 
corresponding to a non-constant wave action flux. On the other hand, the later regime displays a quasi-steady behavior 
consistent with an inverse cascade with a nearly spectrum. This regime exhibits an almost constant flux of 

wave action towards the large scales (see|3]-b) except precisely near the largest scale of the system (fc « 0). 


IV. SIGNATURE OF A FINITE-TIME SINGULARITY 


The first stage of the evolution appears as the formation in time of a spectrum characterized by a non uniform flux 
of wave action from the short to the long scales, as shown in figure ([3]). This flux fills the spectrum from the large k 
towards small fc, tending to a steady power-law spectrum with an almost constant flux of wave action Q (Fig. [5]-b). 

It is important to notice that this built in time spectrum is of finite capacity that is UkcPk < oo (taking 
Uk oc with a ^ 4/3. Therefore one expects, assuming a constant injection of wave action in the injecting domain 

around the formation of such spectrum in finite time. This situation is in fact similar to the self-similar formation 
of a condensate of weakly classical nonlinear waves [ll|) MMM and we shall characterize quantitatively this self¬ 
similar dynamics. To do that, we compute the characteristic length scale involved in the self-similar process, via the 
negative moments (typically n < —2 later on) of the spectral distribution [I[: 


(k-) 


k'^nk{t)cPk ]/{ nk{t)d?k 


(15) 


This allows to define characteristic wave numbers of the spectrum through Fig. |4]shows these characteristic 

wave number to the power 2/3 computed numerically for different moments from n = —2ton = —Tat short times 















(< < 80000 time units). We observe that the different curves exhibit a linear decrease below a critical time suggesting 
the singular behavior for the critical wavenumber of the spectrum: 

ko{t) ^ (t* — with Ri 65000. (16) 

It is the signature of a finite time singularity that would be present if the asymptotic spectrum would be filled with 
a constant wave action flux. 



FIG. 4: The evolution of the characteristic wave number [ko{t)]^^^ computed through the n-th moments of the distribution 
(nsj, for consecutive n ranging from n = —7 up to n = —2 (labeled explicitly on the figure). The simulation conditions are as 
for other figures. The straight lines correspond to a linear fit Kn{t» — i). Notice that almost all moments vanish near an unique 
critical time. The corresponding values for this time are: t* = 85971.8 for n = —7, t, = 76926.9 for n = —6, t, = 69154.8 for 
n = —5, t* = 64971.5 for n = —4, t* = 64173.4 for n = —3, t, = 65376.9 for n = —2. Despite the inaccuracy of the higher 
order moments (-7,-6) all other critical times are around t* « 65000, indicating the independence of t* with order n (Notice 
that the range of the temporal axis is different from the one of Figs. [2]and[3). 


This singular behavior suggests a self-similar solution of the form [Hill: 

- '0 ■ 47 ) 

From relation (1161) . we obtain /3 = 3/2. The parameter a is settled assuming that wave turbulence theory is valid, 
so that the self-similar solution (13 should obey the kinetic equation (|8]). a = 2 is then the only possible choice to 
balance the l.h.s. and the r.h.s terms in the kinetic equation ([5|). Finally, the function satisfies an autonomous 
equation, which reads: 


^(/(s, r) = ^2(/(s, r) -h t)] (18) 

where s = fc(t* — is the self similar variable, r = log(t* — t), and C is formally the same collisional operator of 

dH), but with the scaling variable s instead of k. 

Thus the self-similar function (j){s,T) follows an integro-differential equation (I18|) . with the boundary condition at 
the origin, (/)(0,t) = 0, and with the asymptotic behavior 

= ( 19 ) 

for s —>■ oo and r —> —oo {t —>• t*). The condition \ = 2 — Zv ensures that, in this limit, the tail of the spectrum does 
not depend on time, as is observed in Fig. [S] Notice that one can always re-scale of the nonlinear equation (1181) to 
settle the pre factor in (fTOl) to unity. 

Equation (1181) represents in fact a nonlinear eigenvalue problem for v, which indicates the power law of the spectrum 
at large wavenumber. Such problems are difficult to solve analytically and even numerically since no systematic 
approaches exist PH- Here we will develop an indirect method providing an approximate value only for v. 

Using the relation between the theoretical values of a = 2 and /3 = 3/2, we can rescale the spectra at different times 
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following the self-similar formula (II3 by plotting as a function of s = k/ko{t) taking ko{t) = (k 

A good collapse of the spectra into a single universal curve is then observed in Fig. [5] particularly at large wavenumber 
s. In fact, the envelope of all the rescaled curves defines the function 



FIG. 5: Plot (log-log scale) of the rescaled spectra as function of the rescaled variable s = k/ko{t), with fco(t) = 

This choice shows a very good collapse of all curves in the large k limit. The spectra plotted correspond to times 
ranging from t = 40000 (colored in red) up to t = 60000 (colored in violet) units. 


4>oo{s) = lim 4>{s,r). 

r—¥ — oc 

Averaging the curves near t, for the same simulations but with different system sizes L = 256, L = 512 and 
L = 1024, we obtain a single curve with better resolution for 0oo(s), as shown on Fig. [Hll-a). Then, seeking the 
exponent i/ such that s^'^cjioois) —f 1 for large s, the best fit gives v k, 0.873, which is significantly different (higher) 
than the theoretical value z/ = 2/3 of the inverse cascade eq. (TT^ . Let us emphasize that it is in fact consistent with 
such unsteady regime which fills the spectrum from small to large scales. 

The particular shape of the universal function tpoois) requires a few comments. Firstly, the spectrum decreases near 
s = 0, in agreement with the boundary condition (/(0,t) = 0. Second, as expected, in the self similar variables the 
forcing position in the spectrum tends to s —?► oo, as one approaches the singularity, therefore the forcing only acts 
as a boundary condition in the ultraviolet regime. Finally, the matching region between the inner and outer behavior 
corresponds to the maximum of the function. 


V. VALIDITY OF WAVE TURBULENCE AND THE LATE STAGE REGIME 

The self-similar behavior (EZl), discussed in the previous section, predicts that wave turbulence assumptions will 
not be valid near the finite singularity. In general, wave turbulence theory is no longer valid either because high 
amplitudes of the spectrum are reached at large scale or because of the discrete dynamics of the modes corresponding 
to wave lengths close to the size of the computational domain. 

The nonlinear transition due to high amplitude will appear at small k when, for some wave numbers, the nonlinear 
time scale deduced from equation (|S]) is of the same order as the period of the linear wave • In the present case of 
(HZl) one has that the nonlinear and the linear frequencies scale respectively as: 

oJNhik) ~- — ^ - - 7 and Wk ~ (t* - t)^, (20) 

Uk dt tuf — t 

near t*. Therefore one expects the nonlinearities to be large as f —so that WWT cannot applied anymore. 

In the following we compute numerically two distinct criteria to quantify the ratio between nonlinear and linear 
contributions. First, we computed the ratio between the nonlinear energy, E4 = — J (i(Ax)^ -|- C}) dx dy , and 

the linear energy £’2 = / (|(i9tC)^ + |(^C)^) dxdy , already discussed in @. This is a global criteria which depends 
only on time and indicates the relative importance of both energies in the dynamics. Fig. EK) shows this ratio as 
a function of time. It is observed that E4/E2 is at most of the order of 6 x 10“^, indicating that the non-linear 
contributions (the stretching contributions) to the energy are really small. Incidentally, the maximum of E4/E2 arises 
for t Ri £ confirming the existence of a precursor to a singularity. This can be understood by the following scaling 
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FIG. 6: a) Self-similar universal function (^oo(s) as a function of the self similar variable, s, in linear scale. The data comes 
from the same conditions as Fig. [S] but for three distinct system sizes: L = 256, green dots; L = 512, blue dots; and L = 1024, 
red dots, b) The same self-similar universal function, but in log-log scale. Fitting the exponent u such that cj)oo{s) —>■ 1 for 
large s provides v « 0.873 which is slightly larger than 2/3. 


argument: near the singularity, the quadratic energy would scale like E2 ^ J LOknkCpk ^ (t* — while the fourth 

order energy (the stretching) would follow E 4 ^ f (fk ~ l/(t* — t), hence E 4 /E 2 ^ (t* — t)~^. 

Although this criterium suggests that the nonlinear behavior is globally weak, one cannot ensure that the nonlin¬ 
earities are uniformly weak and in particular that the nonlinearities are effectively small for all scales. A local (in k) 
criteria can be used numerically via the ratio between the linear time scale and the first order nonlinear correction to 
the frequency 0, as stated in eq. (EOl). 

Fig. [7)o) plots the ratio ijj'^^\k)/ijJk, from ([7]), as a function of the wave number at different times. One notices that 
the infra-red behavior of the quotient ui^^\k)/u}k increases significantly when approaching the singularity. Neverthe¬ 
less, it is always less than 10 “^, which implies that the weak amplitude expansion is presumably uniformly valid, even 
near the singularity signature. 

Though wave turbulence theory predicts breakdown of the theory (001) . because ojNL{k)/ojk ^ {t* — one 
observes that direct numerical simulations on the Foppl-von Karman equations (irm do not allow the nonlinearities 
to be of the order of unity. Nevertheless, in the case of a large forcing, such effects have been seen in our numerics 
and it could be then a reason for the breakdown of wave turbulence. 

In conclusion, in the limit of small forcing investigated here, the system does not create strong nonlinearities 
although the singular behavior is cured near t*. This effect comes from the discrete properties of the system which 
becomes relevant at this stage, in particular, to the discrete dynamics of the first modes of the pla te ( the lowest in 
term of frequency) which have to be considered in a modified picture of wave turbulence theory |17l . l33l| . 


VI. DISCUSSION 

Our numerical study reveals that a wave action inverse cascade is built in time and eventually reaches the infrared 
region in finite time through a clearly identifiable self similar process. In this early time regime, the wave system is 
driven by the WWT kinetic equation (|8]), and the dynamics is characterized by a self similar evolution which should 
eventually blow-up in finite time. However, near the singularity, the dynamics is smoothed and the kinetic equation 
is no longer valid: in the small forcing cases, investigated here, the system is governed by the discrete dynamics 
of the largest modes coupled with the continuous spectrum (a discrete breakdown). On the other hand, for larger 
forcing (not studied here), a regularization of the dynamics through the nonlinear breakdown of the WWT is expected 
(nonlinear breakdown). 

For other systems where the wave action is a conserved quantity (for instance for the nonlinear Schrddingier 
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FIG. 7: a) The ratio EilE 2 vs. t, under the same conditions as in Fig. [3]but for three distinct system sizes L = 256, 512 and 
1024 with 512^, 1024^ and 2048^ modes respectively, b) The ratio /ijJh as a function of the wave nnmber k for different 

times steps in log-log scale. The range of time plotted corresponds exactly to the one used in Fig. [5] that is from t = 40000 
time units (colored in red) np to t = 60000 (colored in violet). Both qnantities show the weakness of the nonlinearities in all 
wave numbers and at all times, justifying the validity of the weak turbulence theory. 


equation), a condensate at (or around) fc = 0 forms, changing the post blow-up dynamics Such effects are not 

possible here since the wave action is not conserved by the dynamics. The mode fc = 0 is neutral and remains null with 
time. Nevertheless, as it has been shown for the non-linear Schrodinger equation [s^, the first modes of the systems 
can exhibit an autonomous dynamics. Figure [8]-a) shows precisely the evolution with time of the amplitude for the 
lowest mode of the plates: after the blow-up time t*, the amplitude of this mode grows more or less linearly in time, 
but exhibiting also important oscillations. Finally, for these large times, a stationary regime is eventually reached 
where the spectrum behaves approximatively like Uk for the low wave numbers (fc < 0.2 in Fig. |8]-b) and 

rife ~ \/k^^o.873 £gj. smaller wavelength (0.5 < fc < fc^ in Fig. [5]-b). This stationary regime is not surprising, since 
we expect that the long time behavior does not change the power law built by the self similar evolution (Section IIIII) , 
and long wave modulations transfer wave action toward fc = 0 with precisely the Kolmogorov-Zakharov spectrum, 
Tik ^ which corresponds to the constant wave action flux solution found for C 22 [nfc]. However, notice that 

this regime is not formally steady for the smallest wave numbers, where the dynamics appears rather more like a 
relaxation towards a full stationary dynamics. 




FIG. 8: a) Evolution of the fraction of the first mode. This plot corresponds exactly to the numerical simulation of Fig. [T]and 
the times of the different snapshots are indicated, b) The long-time evolution of the spectra. The time considered are from 
t = 80000 units (colored by red) up to t = 240000 (colored by violet). The spectra only vary for the small wavenumbers where 
it relaxes towards the stationary power law spectrum rife oc k~‘^'' with v ~ 2/3. The line (i) correspond to a 4/3-power law, 
while the line {it) corresponds to a 2 x 0.873 power law which represents better the behavior in the large fc limit. 
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In conclusion, although an elastic vibrating plate does not formally posses a wave action conservation law, an 
undoubtedly inverse cascade of wave action is observed, exhibiting a complex time dependent dynamics. The process 
of formation of such an inverse cascade is ruled by a self similar evolution of the spectra which transfers wave action 
from short wavelength scales up to long wave length scales. Formally the observed self similar solution blows-up leading 
to a singular behavior which is cured in the original system, probably, because of the finite size of the system and the 
role played by the discreteness of the modes. The late evolution of the system is also governed by wave turbulence 
theory, although the discrete nature of the lowest modes modifies the overall picture. This scenario is consistent with 
the formation of a coherent structure which is characterized by the largest modes of oscillation (see Fig. [7^) plus 
small fluctuations. Remarkably, the nonlinear fraction of the energy indicates that this coherent structure makes the 
stretching very small. 
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